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Application  of  Complex  Ray  Theory  in  EM  Scattering  and  RCS  Analysis 


Ruan.  Yingzheng  (University  of  Electronic  Science  and  Technology  of  China. 

Chengdu) 

(Manuscript  received  in  November,  1987.  and  revised  in  October,  1988.  This 
topic  was  sponsored  by  the  National  Educational  Foundation.) 

Ar^-rv^cr :  ^  The  complex  ray  theory  provides  a  simple  and  useful  approach 

for  the  analysis  of  the  propagation  and  scattering  of  HF  EM  waves  in 

various  complicated  environments.  This  paper  deals  with  complex  ray  fields 

in  free  space,  reflection  and  refraction  at  a  single  interface,  multiple 

reflections  and  transmission  through  a  layered  medium,  total  reflection  and 

critical  reflection  at  a  planar  interface,  diffraction  at  an  edge,  complex 

ray  expansion  and  superposition,  as  well  as  RCS  calculation  by  complex  ray 

method.  In  addition,  some  numerical  results  are  given.  ■')* 

A  •'  '•*  *  ‘  /Stsf*  /*  + 
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I.  Complex  Ray  in  Free  Space  (1~'5) 

Complex  ray  theory  is  a  theory  developed  in  the  70 ’s  to  obtain  an 
approximation  solution  to  the  high  frequency  wave  field  problems.  In  the 
80' s.  due  to  the  rapid  progress  in  the  computer  technology,  the  theory  of 
complex  ray  theory  has  received  wide  application  in  the  electromagnetic- 
wave  related  fields  such  as  laser,  optical  fiber,  antenna,  and  microwave 
and  the  elastic  wave  fields  such  as  water  sound  and  earthquake  wave.  The 
purpose  of  this  paper  is  to  discuss  the  application  of  complex  ray  theory 
in  the  study  of  high  frequency  electromagnetic  scattering  and  analysis  of 
radar  cross  section. 

According  to  analytical  expansion  principle  of  complex  function,  if  the 
magnitude  of  the  imaginary  part  of  the  coordinate  of  the  waver  soui ce  is 
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given,  the  complex  source  point  with  complex  coordinate  can  be  obtained. 
These  rays  propagating  in  the  complex  space  from  the  complex  source  and  are 
called  the  complex  rays.  Following  the  traditional  ray  theory  and  methods, 
the  trace  of  the  complex  rays  can  be  followed  and  the  ray  field  can  be 
calculated. 

In  order  to  simplify  the  problem,  the  discussion  in  this  paper  was  based 
on  a  two-dimensional  structure  Figure  1  shews  the  free  space.  A 

wave  vector  b  was  defined  at  the  source  of  the  polarized  linear  electrical 
current  S(ys,zs)  along  the  x  direction.  The  parameter  b  was  called  the  wave 
width  parameter  and  the  directional  angle  a^  was  called  the  wave 
directional  angle  of  the  vector  b  (see  figure  1).  The  coordinate  of  the 
complex  source  was  defined  as  s(ys.  zs) 


Using  the  cylindrical  Bogarin  function  far  expansion,  the  field  expression 
at  the  observation  point  P(y,z)  of  the  two-dimensional  complex  source 
(ignoring  the  time  factor  exp(-jwt))  is 


f  «  ~ 
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where  the  complex  distance  is 


In  the  real  space,  the  plane  that  pa— a  through  S  point  and  is 
perpendicular  to  the  vector  b  is  defined  as  the  caliber  plane  of  the 
complex  source  (see  figure  1) .  On  this  caliber  plane,  the  oceplex  source 
point  field  has  the  uniform  phase  distribution  and  the  Gauss i an-type 
amplitude  distribution  and  its  remote  radiation  field  is  also  a  Gaussian- 
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Figure  i:  Two-dimensional  complex  source  point  and  its  caliber  field 
distribution 


key:  1  -  caliber  circle 


type  fixed-directional  radiation  wave.  Therefore,  it  is  appropriate  to  use 
the  complex  source  point  field  to  simulate  the  concentrated  waves  such  as 
the  mono-pulse  laser  wave  or  the  radar  antenna  major  wave  front.  Based  on 
the  conclusion  of  the  real-space  ray  optics  and  after  simple  analysis  and 
expansion,  the  transmission  and  scattering  of  the  wave  field  can  be  solved. 
Since  the  physical  optical  integration  or  the  plane  wave  spectrum 
integration  of  the  caliber  field  is  avoided,  the  numerical  process  was  made 
even  simpler. 

In  this  paper,  we  will  briefly  introduce  the  various  methods  of  the 
complex  ray  theory,  including  the  complex  ray  tracing  method,  the  near-axis 
approximation  method,  the  agglomeration  ray  method,  the  complex  ray 
diffraction  method,  the  complex  ray  expansion  method,  and  the  application 
of  these  methods  in  the  calculation  of  the  scattering  field  and  analysis  of 
the  radar  cross  section. 

II.  Reflection  and  Refraction  of  Complex  Rays [6-7] 

Considering  the  two-dimensional  randcxn  curved  surface  f(y,z)K)  shown  in 
ficrure  2.  for  the  complex  source  S  defined  in  equation  (1)  and  the 
observation  Doint  P.  a  complex  ray  which  passes  through  a  point  A  on  the 
curved  surface  and  reaches  point  P  can  be  obtained  through  iteration 
process.  This  ray  should  satisfy  the  complex  space  Feima  principle  with  the 
complex  light  distance  as  the  limit  value 

ste&m-., .  ,4, 

where  Lg  and  L2  are  the  complex  distances  for  the  incident  and  the 
reflected  paths  and  a  is  the  complex  reflected  angle  (see  figure  2). 
Therefore,  the  field  at  point  P  can  be  written  as 
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where  R(@f)  and  6(0)  are  the  complex  reflection  and  complex  diffusion 
functions.  These  functions  can  be  obtained  by  substituting  the  incident 
angle. 0  in  the  real  space  with  the  complex  incident  angle  8.  The  method 
which  uses  the  iteration  technique  to  trace  the  complex  ray  parameter  (such 
as  the  complex  reflected  angle  CL)  is  called  the  complex  ray  tracing  method. 
Obviously,  tremendous  computer  time  is  required  in  the  tracing  process. 

In  order  to  simplify  the  numerical  process,  the  near-axis  approximation 
can  be  used. [6,7]  it  only  requires  the  tracing  of  the  real  trace  along  the 
wave  axial  line,  then  the  complex  ray  field  Eq  at  the  axial  observation 
point  P0  can  be  calculated  approximately  and  then  the  complex  ray  field  at 
the  near-axis  observation  Doint  P  can  be  obtained  through  the  micro¬ 
disturbance  correction  of  complex  phase  and  complex  amplitude 


where 
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(6) 
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The  %  in  equation  (6)  denotes  the  complex  phase  correction  and  vp  and  vp 
are  the  corrections  for  the  reflection  and  diffusion  functions.  They  are 
all  functions  of  the  near-axis  distance  d  and  can  be  expanded  by  Taylor’s 
series.  (In  Taylor's  series  expansion,  only  the  first  term  should  be 
considered  for  approximation  calculation.  [6,7}) 

For  the  complex  rays  that  penetrate  through  the  interface  and  into 
another  medium,  the  complex  expansion  should  be  based  on  the  Snare's  law 
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Ficrur e  3 


Near-ax is  approximation  of  coup lex  rays 
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and  the  results  of  the  refraction  of  complex  rays  should  be  similar  to 
those  of  the  reflection  of  complex  ravs. 

III.  Continuous  Refraction  and  Multiple  Reflection  of  Complex  Rays^®-^ 

In  the  layered  homoaeneous  medium  shown  in  figure  4.  the  complex  rays 
can  reach  the  observation  point  after  multiple  reflections  and  refractions. 
Therefore .  the  total  field  at  point  P  should  be  the  summation  of  all  the 
possible  complex  ray  fields 


f. 


(8) 


where  denotes  the  complex  ray  field  which  reaches  point  P  after  M  times 

of  reflections  between  the  two  layers  (see  figure  4).  The  calculation  can 
be  proceeded  with  either  complex  ray  tracing  method  or  the  near-axis 
approximat ion  method f-' 


(9) 


wnere  T-t ; f 2  is  the  coefficient  of  refraction  of  this  double— refraction. 

.  **' 

*?„»■  nr.  is  the  total  coefficient  of  reflection  after  M  times  of 

•  -  1 

reflections.  5w***nl  2.  is  the  total  diffusion  coefficient,  and  rf*  Lm+  l, 

*- I 

i«  the  total  complex  light  distance.  They  can  all  be  solved  from  the 
complex  ray  paths. [8-10] 

In  order  to  simplify  the  summation  process  of  the  infinite  series  in 
equation  (8) ,  the  residual  terms  from  the  M-th  term  can  be  treated 
together.  If  the  thickness  of  these  layers  is  uniform  or  if  the  incident 
ray  is  near  the  normal  vector  of  the  curved  surface,  then  the  residual 
terms  can  be  expressed,  after  Fourier  transformation  and  integration 
approximation  calculation,  ast-^ 
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Comparison  of  equations  (8)  and  (10)  shows  that  the  summation  of  these 
residual  terms  cam  be  regarded  as  the  contribution  from  the  M-th  complex 
ray.  It  would  only  require  to  replace  the  transmission  coefficient  T  in 
equation  (8)  by  the  treinsmission  coefficient  of  the  agglomeration  of  the 
rays  Tc.  Therefore,  equation  (8)  can  be  re-written  as  the  mixed-field 
structure  of  the  ordinary  rays  and  the  agglomeration  rays 


-  1  ~  K 

£=  r  Eh+Ec 

M-l 


(11) 


The  numerical  results  show  that  a  satisfactory  scattering  field  result, 

from  a  engineering  standpoint,  can  be  obtained  with  M-0  (no  ordinary  rays, 
only  the  agglomeration  rays)  or  M=1  (one  ordinary  ray,  and  the 
agglomeration  rays) . 


IV.  Total  Reflection  and  Critical  Reflection  of  Complex  Rays  112-14] 

When  the  complex  ray  field  with  Gaussian  wave  front  is  incident  from  a 
optically  less  dense  medium  to  an  optically  more  dense  medium  and  if  on  the 
interface  (see  figure  5,  c,>e,),  the  incident  angle  0  is  less  than  the  total 
reflection  critical  angle  e,  •  •ur*,v'e,/e1  than  there  exists  only  the 
geometrical  optical  reflection  field  Eg,  irtrtch  is  the  contribution  from  the 
complex  ray  field  to  the  independent  saddle  point  and  is  expressed  in  terms 
of  the  integration  of  the  plane  wave  spectrum  Eg 

%kJ7*th''"9'  t m 

where 

-v'o- >.)*♦  (»♦.,)*  ,»•<*,)><>  <»»> 
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Furore  5:  Total  reflection  and  lateral  waves  of  complex  rays 
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When  the  incident  angle  is  greater  them  the  critical  angle,  the  total 
reflection  field  of  the  complex  ray  should  include  both  the  contributions 
to  the  independent  saddle  point  Eg  and  to  the  independent  supporting  point 
E^] 12-14]  from  the  geometrical  optical  field  EJ-f  and  the  complex  lateral 
wave  field  £) :  namely 


(14) 


where 


£,=£.«- 


(15) 


where  L^.  To-  and  L3  denote  the  transmission  paths  of  the  lateral  waves  as 
shown  in  figure  5.  Obviously,  when  0-X)c.  I^-X)  and  equation  (15)  gives  the 
quantitative  results  of  the  spreading  field  and  there  exists  a  transition 
region  of  the  complex  ray  reflection  field.  In  order  to  make  the 
singularity  of  the  transition  field  smooth,  a  uniform  equation  which 
includes  the  transitional  function  should  be  used.  In  other  words,  the 
asvmptotic  contribution  when  the  saddle  point  is  close  to  the  supporting 
point  should  be  expressed  as [12-14] 


E-El*E  i  (16) 

where  E'a  is  given  by  equation  (12)  with  R(6)«l,  is  the  field  of  the 
transitional  region 


e  f  *  i*  ,  t r  ft*h/tT} 

•iaeccMie,i- 


(17) 
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where 


■ip(flf  -  ft) 
•in2  6  c 


(18) 


r  = 


I  /  c — ?  +  »in? 

g\  e©«5Sc 


(19) 


(20) 

where  Di/2(x)  Is  the  1/2  order  parabolic  column  function  of  the  complex 
parameter  x.  The  ranere  of  the  transitional  region  is  determined  by  the 
magnitude  of  the  complex  parameter.  Therefore,  the  parameter  Sc  can  be 
defined  as  the  "numerical  distance"  of  x  and  the  numerical  results  show 
that  x=2  is  a  criterion  for  the  judgement : [12,13]  when  x>2  then  the  general 
equation  (14)  can  be  used  and  when  x< 2  then  the  uniform  equation  (16) 
should  be  used. 

/ 

Besides  the  total  reflection  transitional  region,  singularities  of  the 
complex  ray  field  also  exist  in  the  reflection  interface,  the  shadow 
interface,  the  focus  line,  and  the  focus  point  regions.  When  this  is  the 
case,  an  appropriate  special  function  should  be  used  for  the  smoothing- 
calibration  of  the  transitional  function  propagating  field  as  described 
above. 


V.  Diffraction  of  Complex  Rays^J 

Diffraction  of  complex  rays  occurs  at  the  edge,  tip,  or  smooth  convex 
surfaces  of  the  structure.  Figure  6  shows  a  Gaussian  wave  incident  upon  art 
ideally  conducting  half-plane  structure.  If  the  complex  source  point  field 
is  used  to  express  the  Gaussian  wave  and  if  the  extension  of  the 
aoplicability  of  the  theoretical  equations  of  geometrical  diffraction  of 
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Figure  6:  Blue  diffraction  of  complex  rays 
key:  1  -  metallic  half  plane 
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linear  source  to  the  complex  space  is  made,  then  the  complex  ray- 
diffraction  field  can  be  obtained. 


MM 


(21) 


where  0O  is  the  complex  incident  angle  at  the  edge  of  the  half  plane,  0  is 
the  diffraction  angle  at  the  observation  point,  and"ft(9,0o)  is  the  complex 
diffraction  coefficient 


5<9j')m-7jzzir 

(22) 

Obviously,  if  the  observation  point  is  located  within  the  transitional 
region  of  the  reflection  interface  or  the  shadow  interface  (§^->0) ,  then 
the  above  equation  loses  its  applicability.  When  this  is  the  case,  then  in 
the  wave  spectrum  integration  formula,  the  saddle  point  is  very  close  to 
the  supporting  point  and  the  simple  geometrical  diffraction  theory  can  not 
be  applied  and  the  uniform  diffraction  theory  and  related  equations  should 
be  used  for  the  analysis  of  the  complex  space. 

The  complex  ray  theory  has  been  used  in  the  analysis  of  various 
diffraction  problems  involving  waves  such  as  the  curved-edge  diffraction  of 
the  antenna  of  ref  lectors  .c's:i 


VI.  Expansion  of  Complex  Rays  and  Superposition  of  Complex  Source 
Points  (16-1.9) 

Based  on  the  principles  of  complex  ray  theory,  we  can  extend  Idle  •' 
Wilkings-Feiner 's  principle  and  its  Gilhoff  mathematical  expression  to  the 
complex  space[16,17j  and  develop  a  new  wave  expansion  method. 
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Considering  the  two-dimensional  structure  shown  in  figure  7  and  assuming 
2-2’  denotes  the  Wilkings  plane,  then  the  wave  source  distribution  function 
on  this  plane  f(y')  is 


_  oo<y'<oo  (23) 

The  field  at  observation  point  P  can  be  expressed  by  the  Gilhoff  equation 


if  the  integration  expression  of  the  source  distribution  is  known 


(24) 

where 


f>  “  (*-*')* 


Obviously,  f(y')  is  the  weighted  function  of  the  source  distribution.  After 
Fourier  transformation,  equation  (24)  can  be  expressed  in  terms  of  the 


integration  of  plane  wave  spectrum 


£*if“  — ■?—), 

2J-V  i-T*  (26) 

where  F(q)  is  the  plane  wave  spectrum  function  corresponding  to  the  field 
distribution  f (y1 ) . 

If  the  Wilkings  plane  is  extended  to  the  complex  space  or  if  the  sub- 
wave  source  is  extended  to  the  complex  Wilkings  source,  then  from  (24)  one 
can  obtain 

W\" <**)*</ >4?  i 

•  -—w-r - r-—  (27) 

where  p  is  given  by  equation  (3).  W(y')  denotes  the  weighted  function  of 
the  superposition  of  complex  source  fields  (sometimes  called  the  expansion 
coefficient  of  complex  ray) .  Since  complex  source  field  is  a  Gaussian-type 
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r igure 
kev :  1 


% 


/ :  Expansion  of  complex  rav 
-  Walkings  plane 
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wave  field  with  fixed  direction  and  the  energy  is  concentrated  around  the 
axis  of  the  wave  front,  the  integration  limits  in  equation  (27)  can  be 
reduced  from  (-00.  00)  to  a  smaller  range  Also  the  discrete 

complex  source  field  can  be  used  to  replace  the  continuous  integration  in 
equation  (27)  and,  therefore,  simplify  the  numerical  process.  Furthermore, 
from  the  complex  expansion  of  the  sub-wave  sources,  the  characteristics  of 
the  field  transitional  region  can  be  changed,  112,13]  the  solution  of 
the  transitional  field  region  can  be  accomplished  through  the  choice  of 
appropriate  parameters  and  the  direct  application  of  the  superposition  of 
complex  source  fields  and  the  singularities  of  the  field  can  be  made  smooth 
without  the  introduction  of  complicated  transitional  functions. 

In  order  to  apply  the  complex  ray  expansion  equation  (27) ,  the  weighted 
function  W(y' )  must  first  be  solved  and  then  the  appropriate  free  parameter 
for  expansion  should  be  chosen.  These  parameters  are  the  ones  such  as 
Wilkings  plane  position  z' .  modulus  and  direction  of  wave  vector  b, 
expansion  range  (yi'.y2')«  discrete  source  distance  hy’,  and  computer 
threshold  value  a,-, .  There  have  been  detailed  discussion  on  the  complex  ray 
expansion  of  cylindrical  waves, and  the  results  have  been  applied  to 
the  reflection  at  plane  interfaces. 

When  b-X),  the  complex  source  field  transforms  to  the  real  source  field, 
and  the  expansion  of  complex  rays  transforms  to  the  general  Wilkings 
principle,  when  b->oo,  the  complex  source  field  transforms  to  the  plane 
wave  radiation  and  the  integration  of  complex  source  field  transforas  4&o 
the  integration  of  plane  wave  spectrum.  Hence,  it  can  be  shown  that  tffe 
expansion  of  complex  source  points  is  a  more  general  farm  of  womb  field 
expansion  than  the  source  distribution  integration  and  plane  wmve  spectrum 
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integration.  The  latter  two  cases  can  be  viewed  as  the  special  cases  of  the 
first  case. 


7.  Complex  Ray  Analysis  of  Radar  Cross  Section [20-22] 

Based  on  the  theoretical  definition  of  target  radar  cross  section  (RCS) 


a  « 


1  i*  tslf 


(28) 


one  can  see  that  the  calculation  of  the  radar  cross  section  can  be 
regarded  as  the  calculation  of  the  scattering  wave  amplitude  at  the  receipt 
point  lEgl  if  the  incident  wave  amplitude  lE^I  is  given.  Therefore,  in 
principle,  all  the  methods  used  to  analyze  scattering  of  electromagnetic 
field  can  be  used  in  the  calculation  of  radar  cross  section.  Since  the 
Gaussian  wave  of  the  complex  source  is  a  very  good  simulation  of  the  major 
wave  packet  of  the  radar  antenna,  the  complex  ray  theory  can  be  used  to 
obtain  the  scattering  field  at  the  target  and  the  single-station  or  double- 
station  radar  cross  sections. 

Figure  8  shows  the  calculational  results  of  the  cross  section  of  an 
intermediate  medium  single-station  radar  cross  section. Except  for  the 
complex  ray  reflection  and  refraction  on  the  surface  of  the  medium, 
diffraction  occurs  at  the  edge  of  the  rounded  plate  and  the  contribution 
from  the  diffraction  should  be  included  in  the  calculation  if  the  diameter 
of  the  rounded  plate  is  not  large.  In  figure  8,  the  solid  line  is  the 
integration  results  based  on  physical  optics  and  dotted  line  is  the 
calculational  results  based  on  complex  ray  theory.  The  distance  of  the 
observation  point  is  taken  as  R-IO12*.  Tran  this  figure,  we  can  see  that 
when  the  incident  angle  is  less  than  50°,  within  the  range  of  the  major 
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Ficrure  8:  The  medium  circular  Dlate  single-station  RCS  calculational 

]  tc 

kev:  5  -  incident  angle 

6  -  Radar  cross  section  (db  square  meter) 

7  -  parameters  used  in  calculation 

8  -  solid  line  —  Dhvsical  optics  method 

9  -  points  —  xxx  complex  ray  method 
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petal  and  the  near  petal,  the  difference  In  the  results  based  on  these  two 
methods  Is  less  than  3db. 

For  the  large-diameter  concave  cavity  structure  such  as  the  gas  inlet 
and  outlet  duct,  if  the  radar  cross  section  is  calculated  based  on  the 
traditional  geometrical  optics  methods,  then  it  is  required  to  form  a  set 
of  incident  rays  on  the  diameter  surface  and  the  distance  between  these 
rays  should  be  less  than  0.2X.  Obviously,  a  lot  of  computer  time  would  be 
necessary  to  analyze  the  reflection  of  a  few  tens  even  hundreds  of 
geometrical  rays.  If  the  analysis  is  conducted  in  the  guided-wave  fashion, 
then  since  the  dimension  of  the  wave  guide  is  far  greater  than  that  of  the 
radar  waves,  the  expansion  and  analysis  of  a  few  hundreds  high  order  modes 
is  still  not  easy.  If  the  complex  ray  theory  is  used  in  expansion  and  near¬ 
axis  approximation  is  used  in  simplification  of  the  problem,  then  tracing 
of  a  few  reflection  paths  of  the  complex  rays  will  be  sufficient  to 
establish  the  field  distribution  on  the  lip  surface  of  gas  inlet  duct  or 
the  field  strength  at  the  reception  point  can  be  determined  and  then  the 
radar  cross  section  of  the  target  can  be  calculated.  The  RG5  analytical 
results  of  a  curved  gas  inlet  duct  has  been  completed.  The  results 
based  on  complex  ray  expansion  were  compared  with  those  based  on  other 
methods. 

The  radar  cross  section  of  antenna  and  antenna  cap  can  also  be 
calculated  by  the  complex  ray  theory. ^2]  The  structures  of  antenna  and 
antenna  cap  with  low  RCS  can  be  optimized  in  this  way. 

VIII.  Concluding  Remarks 

The  complex  ray  theory  has  been  proved  to  be  an  efficient  way  of 
analysis  of  high  frequency  field  scattering.  This  paper  discusses  the 
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application  of  this  theory  under  the  conditions  of  uniform  and  layered 
medium.  These  are  just  some  simple  conditions.  However,  its  application  is 
not  limited  by  these  conditions.  In  the  non-uniform,  dissipated,  or 
directional ly-conducting  media,  the  complex  ray  theory  still  has  very 
important  applications.  Using  the  complex  ray  theory,  not  only  can  the 
radar  cross  sections  of  the  metallic  parts  such  as  the  body  and  wings  of 
the  vehicle  be  calculated,  the  reflection  and  transmission  of  the  radar  cap 
and  cabin  hoods  can  also  be  analyzed.  Not  only  can  the  scattering 
contribution  due  to  the  key  scattering  parts  such  as  antenna  and  gas  inlet 
duct  be  calculated,  but  also  the  effect  of  RCS  reduction  after  application 
of  wave-absorbing  materials  can  be  analyzed.  So  far,  this  topic  has  been 
under  intensive  study  worldwide  and.  in  some  areas,  significant 
accomplishments  have  been  made. 
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